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Abstract 



A microscopic angular momentum projection after variation is used to describe 
quadrupole collectivity in 30 ' 32 ' u Mg and 32,34,36,38 rp^g jjartree-Fock-Bogoliubov 
states obtained in the quadrupole constrained mean field approach are taken as in- 
trinsic states for the projection. Excitation energies of the first 2 + states and the 
B(E2,0 + — » 2 + ) transition probabilities are given. A reasonable agreement with 
O I available experimental data is obtained. It is also shown that the mean field picture 

of those nuclei is strongly modified by the projection. 
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Neutron-rich nuclei with A" « 20 are spectacular examples of shape coexistence 
between spherical and deformed states. Experimental evidence for an island 
of deformed nuclei near A^ = 20 has been found in the fact that 31 Na and 
32 Na are more tightly bound than could be explained with spherical shapes 
[1]. Additional support comes from the unusually low excitation energy of 
the 2 + state in 32 Mg [2]. A large ground state deformation has also been 
inferred from intermediate energy Coulomb excitation studies [3] in 32 Mg. 
Quadrupole collectivity of 32-38 Si has also been studied in [4]. Very recently 
this region has been the subject of detailed experimental spectroscopic studies 
at ISOLDE [5]. From a theoretical point of view, deformed ground states 
have been predicted for nuclei with N ^ 20 [6-8]. In those calculations the 
rotational energy correction is the essential ingredient for the stabilization 
of the deformed configuration. On the other hand, some calculations have 
predicted [9-12] a spherical ground state for 32 Mg but it has also been found 
[12] that deformation effects may appear as a result of dynamical correlations. 
Some shell model calculations, even with restricted configuration spaces, have 
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been able to explain the increased quadrupole collectivity at N = 20 as a 
result of neutron 2p — 2h excitations into the fp shell, see for example [13,14]. 
Recently, a mean field study has explored the suitability of several Skyrme 
parameterizations [15] in the description of this and other regions of shape 
coexistence. 

The mean field description of nuclei is usually a good starting point as it 
provides a qualitative, and in many cases quantitative, understanding of the 
nuclear properties. This is the case when the mean field solution corresponds 
to a well defined minimum. However, in regions of shape coexistence where 
two minima are found at a comparable energy, the correlation effects stem- 
ming from the restoration of broken symmetries and/or collective motion can 
dramatically alter the energy landscape thus changing the mean field picture. 
For this reason, we have included in our mean field calculations the effects 
related to the restoration of the broken rotational symmetry by performing, 
for the first time with the Gogny force, angular momentum projected calcu- 
lations of the energies and other relevant quantities. The reason for choosing 
to restore rotational symmetry is that the zero point energy associated with 
this restoration is somehow proportional to deformation and ranges, in this 
region, from a few KeV for nearly spherical configurations to several MeV for 
well deformed ones. This energy range is comparable to the energy differences 
found between different shapes in nuclei of this region. Therefore, in addition 
to the mean field results, both angular momentum projected 1 = and 1 = 2 
surfaces were computed for the nuclei 30 > 32 > 34 M<7 and s 2 ' 34 . 36 * 38 ^ and angular 
momentum projected transition probabilities B(E2, + — > 2 + ) among different 
configurations. 

The calculation proceeds in two steps: in the first one we perform a set of 
constrained Hartree- Fock- Bogoliubov (HFB) calculations using the D1S pa- 
rameterization [16] of the Gogny force [17] and the mass quadrupole operator 
Q20 = z 2 — |(x 2 + y 2 ) as the constraining operator in order to obtain a set 
of "intrinsic" wave functions | 4>(q2o)}- The self-consistent symmetries imposed 
in the calculation were axial symmetry, parity and time reversal. The two 
body kinetic energy correction was fully taken into account in the minimiza- 
tion process. On the other hand, the Coulomb exchange term was replaced by 
the local Slater approximation and neglected in the variational process. The 
Coulomb pairing term as well as the contribution to the pairing field from 
the spin-orbit interaction were neglected. A harmonic oscillator (HO) basis of 
10 major shells was used to expand the quasi-particle operators and the two 
oscillator lengths defining the axially symmetric HO basis were kept equal for 
all the values of the quadrupole moment. The reason for choosing the basis 
this way was that we wanted a basis closed under rotations (i.e. an arbitrary 
rotation of the basis elements always yields wave functions that can be solely 
expressed as linear combinations of the elements of the basis) in order to avoid 
the technical difficulties discussed in [18] when a non-closed basis is used. In 
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the second step we compute the angular momentum projected energy for each 
intrinsic wave function | 0(?2o)) obtaining in this way a set of energy curves 
Eifao) for each value of I = 0, 2, . . . The minima of each curve provide us 
with the energies and wave functions of the / = + , 2 + , . . . yrast and isomeric 
states. 

The theoretical background for angular momentum projection is very well de- 
scribed in [19,20] and therefore we will not dwell on the details here. However, 
a few remarks concerning the peculiarities of our calculation are in order: 
first, and due to the axial symmetry imposed in the HFB wave functions, the 
angular momentum projected energy is given by 

v< \ Jo d(3sen(3dU(3){<P(q w ) \ H'[p {r)]e-^y | 0(g 2O )) 
Jo 2 d/3sen/3d^(P)(cl>(q w ) | e~^y | 0(g 2O )) 



with H'[p p (r)} = H[p (3 (r)]-\ w (N n -Z)-X u (N u -N). The term -A^-Z)- 
\ U (N U — N) is included to account for the fact that the projected wave function 
does not have the right number of particles on the average. The previous term 
would correspond to the application of first order perturbation theory if the 
chemical potentials used were the derivatives of the projected energy with 
respect to the number of particles. In our calculations we have simply used 
the chemical potentials obtained in the HFB theoryQ. This simplification is 
justified by the fact that the deviations induced in the number of particles due 
to the angular momentum projection are always small and so are their effects 
on the projected energies. 

For the computation of the matrix elements of the rotation operator in a HO 
basis we have used the results of ref. [21]. 

Another relevant point to be discussed is the prescription to use for the density 
dependent part of the Gogny force. In the calculation of the energy functional 
E[4>] = (0| H \ 4>) the density appearing in the density dependent part of the 
force is simply p(r) = (0| p \<p) rendering the energy a functional of the density 
and the pairing tensor but with a dependence on the density different from 
the canonical quadratic one of the standard HFB theory. On the other hand, 
the energy overlap E[<fi, (f)'} = (<f)\H \<p') / {4>\4>') can De evaluated using the ex- 
tended Wick theorem. The final expression is the same as the HFB functional 
E[(f>] but replacing the density matrix by pij = (<p\ c+q / (4>\<j)') and the 
pairing tensor by Rij = (<f>\ CjCi / {4>\4>') an d i^ij = (<f>\ CiCj \<p') / {4>\4>')- As a 
consequence, it seems rather natural to use the density p (r) = (cj>\ p \(j>') / ((f)\4>') 



1 This recipe has been previously used [22] in the context of angular momentum 
projection and in Generator Coordinate Method (GCM) calculations [23]. In both 
cases it has been found that the present recipe works very well. 
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Fig. 1. The HFB and projected energies as a function of g20 for the nucleus 34 Mg. 
See text for further details. 

in the evaluation of the density dependent term of the force in E[cj), (/>']. In our 
case, this leads to the introduction of a density dependent term depending 
on p(r,/3) = (0| pe~ tf3Jy \<p) / (<fi\ e~ ll3Jy \<f>). This density dependence seems to 
yield to bizarre consequences like having a non-hermitian and non rotationally 
invariant hamiltonian. These apparent inconsistencies can be overcome if we 
think of a density dependent force, not as an operator to be added to the 
kinetic energy in order to obtain a hamiltonian, but rather as a device to get 
energy functionals like E[(f)\ and E[(j), </>'] with the property of yielding an en- 
ergy that is a real quantity and independent of the orientation of the reference 
frame. The density dependence just mentioned fulfills these two requirements 
as can be readily checked. In addition, when the intrinsic wave function is 
strongly deformed and the Kamlah expansion can be used to obtain an ap- 
proximate expression for the projected energy (the cranking model) the above 
density dependence yields the correct expression for the angular velocity uj 
including the "rearrangement" term [24]. A more elaborated argumentation in 
favor of the density dependence just mentioned will be given elsewhere. 

As an example of the results obtained we show in Fig. 1 the HFB and projected 
energies as a function of ^20 for the nucleus 34 Mg. In contrast with the HFB 
result, the 1 = energy surface shows two pronounced minima in the prolate 
and oblate side which are rather close to each other in energy being the prolate 
minimum slightly deeper than the oblate one. Therefore, it is difficult to assign 
a given character to the 1 = state until a configuration mixing calculation 
is performed, although it is very likely that the predominant configuration 
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for the 1 = state is going to be the prolate one. For 1 = 2 there is a well 
developed prolate minimum. Let us also mention that for configurations with 
a (720 value close to zero (i.e. close to the spherical configuration g 2 o — which 
is a pure 1 = state) it is very difficult to compute the 1 = 2 projected energy 
due to numerical instabilities related to the smallness of (0(?2o)| Pqq |<K<?2o))- 

In the inset of Fig. 1 we have plotted the energy difference E RO t{I) = E HFB — 
Ej as a function of (720 for / = (full line) in order to compare it with the 
rotational energy correction E^q T = {j^j I Jy often used in mean field calcu- 
lations (dashed line). The Yoccoz moment of inertia Jy has been computed, 
as it is usually done, in an approximate way by neglecting the two body quasi- 
particle interaction term of the hamiltonian (the same kind of approximation 
yields to the Inglis-Belyaev moment of inertia instead of the Thouless-Valatin 
one). We notice that E^q T agrees qualitatively well with Erot(0) for (720 val- 
ues greater than 100/m 2 and smaller than — 50/m 2 as expected: these are 
regions of strong deformation where the validity conditions for E#q T to be a 
good approximation to Erot (Kamlah expansion) are satisfied. On the other 
hand, the behavior of E^q T is completely wrong in the inner region. One pre- 
scription to extend the rotational formula to weakly deformed states is the 
one of [25] based on results with the Nilsson model. The prescription multi- 
plies E^q T by a function of (jy} with the property of going to zero (one) for 

Jy^j going to zero (infinity). The resulting rotational energy is also depicted 
in the inset of Fig. 1 (dotted line) and, although the qualitative agreement 
with the exact result improves somewhat, the quantitative one is far from sat- 
isfactory in the nucleus considered. The rotational energy correction formula 
is based on the assumption that the quantity h((3) = ^He~ l13 '^ / (e~ %f5 ^ y ^ can 

be very well approximated by a quadratic function h(/3) ~ h(0) + \h"{0)f3 2 
where h"{0) is related to the exact Yoccoz moment of inertia by the expression 

Jy = — (Jy) /h"(0). It is well known that this assumption is justified for de- 
formed heavy nuclei. However, we have checked that it is not the case for the 
nuclei studied here even for the largest deformations considered. Therefore, 
we conclude that the exact restoration of the rotational symmetry is funda- 
mental for a qualitative and quantitative description of the rotational energies 
in these light nuclei. 

The main outcomes of the calculation are summarized in Fig. 2 where we show, 
on the left hand side panel, the HFB potential energy surfaces for Mg and Si 
isotopes as a function of the mass quadrupole moment. These surfaces have 
been shifted accordingly to fit them in the plot. We observe that only in the 
nuclei 34 Mg and 38 Si we obtain a prolate minimum at deformations of 0.4 
and 0.35 respectively. For the other nuclei, the minimum corresponds to the 
spherical configuration. For all the nuclei considered the energy curves are very 
flat around the corresponding minimum indicating that further correlations 
can substantially modify the energy landscape and therefore the conclusions 
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Fig. 2. HFB (left), angular momentum projected 1 = (middle) and 1 = 2 (right) 
collective potential energy surfaces for the nuclei 30 ' 32 ' 34 Mg and 32 . 34 . 36 . 38 gi. All 
the curves have been shifted by the same amount for each nucleus in all the three 
cases in order to fit them in a single plot. The energy shifts are -36, -30, and -28 
MeV for w^M Mg an( j _23, -13 and -6 MeV for ^ 2 ^e Si respectively. The range of 
quadrupole moments considered roughly correspond to a /?2 range from -0.5 to 0.85. 

obtained from the raw HFB results. 

On the middle and right hand side panels of Fig. 2 we show the angular mo- 
mentum projected 1 = and 1 = 2 potential energy surfaces for all the nuclei 
considered. These surfaces have also been shifted to fit them in the plot. For 
/ = 0, apart from the nucleus 34 Mg that shows a rather clear prolate minimum, 
the general trend for the ground state is to show shape coexistence. For 1 = 2 
we have prolate minima for 32-34 Mg and oblate minima for 32-34 Si whereas 
the other nuclei are examples of shape coexistent structures. The results just 
shown indicate that, for a quantitative description of the ground and 2 + states 
in all these nuclei, a configuration mixing calculation (GCM) using the mass 
quadrupole moment as generating coordinate is needed. In spite of this, we 
present in Table 1 the + — 2 + energy differences for the four possible configu- 
rations with the + in the prolate (P) or oblate (O) minimum and the 2 + also 
in the P or O minimum. The energies written in boldface correspond to the 
predictions obtained by strictly using the criterion of the absolute minimum of 
Ei(q2o) to assign the + and 2 + states. Comparison of these predictions with 
the experimental results indicates a reasonable agreement except for 36 Si. The 
inclusion of configuration mixing will presumably improve the agreement as 
it always yields to a mixed configuration with an energy lower than the en- 
ergies of the states being mixed. Therefore, if the reconfigurations strongly 
mix (shape coexistence) but the 2 + ones do not ( there is a well established 
minima) the + — 2 + energy difference will increase whereas it will decrease 
if the opposite situation takes place. On the other hand, if configuration mix- 
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Table 1 

Excitation energies, in MeV, for the 2 + states. The four possible combinations are 
shown. The numbers in boldface indicate the configuration where the + and the 
2 + correspond to the absolute minimum of the projected energy surface. The exper- 
imental data are taken from [4] (Si) and from [2] (Mg). 





(0+ - 2+) 


( O " 2 p) 


(0+ - 2+) 


(°o " 2 o) 


Exp 


30 Mg 


1.683 


1.681 


2.391 


2.388 




32 Mg 


0.873 


0.235 


4.546 


3.909 


0.885 


3A Mg 


0.753 


-1.206 


5.960 


4.000 




32 Si 


4.475 


5.408 


0.869 


1.803 


1.941 


u Si 


4.088 


4.324 


2.383 


2.619 


3.327 


m Si 


2.845 


2.670 


2.778 


2.603 


1.399 


38 Si 


1.238 


1.481 


1.837 


2.080 


1.084 



ing is important for both the + and 2 + states anything can happen to the 
excitation energy. Therefore, we expect that configuration mixing is going to 
increase the excitation energies in all cases except in 30 Mg and 36,38 Si where 
the behavior is unpredictable. In Table 2 we present the results obtained for 
the B(E2, + — > 2 + ) transition probabilities for the four possible combina- 
tions. As in the previous table, the results obtained by choosing for the + 
and 2 + states the ones corresponding to the absolute minima of the projected 
energies are written in boldface. For the nuclei 32 Mg and u Mg we obtain very 
collective values for the B{E2) which, in the case of 32 Mg, are in rather good 
agreement with the experiment. For both nuclei, we expect a contamination 
of the ground state wave function by the oblate + state that will yield to a 
reduction of the B(E2) values (see column two for the B(E2, 0q — > 2p)) that 
will bring the theoretical predictions in closer agreement with the experimental 
data. For the 32 Si, 34 Si and 38 Si isotopes we underestimate the B(E2) values 
but, presumably, admixtures of the Op — > 2p transition will help to bring the 
theoretical results in closer agreement with the experiment, specially for the 
38 Si nucleus. Concerning 36 Si we can only conclude that a strong Op — > 2p 
component has to be present in the evaluation of the B(E2). 

In Table 3 the HFB and projected ground state energies for the nuclei under 
consideration are shown and compared to the experimental data taken from 
[26]. The inclusion of the zero point energy stemming from the restoration 
of the rotational symmetry clearly improves the theoretical description of the 
binding energies. 

In conclusion, we have computed several properties of neutron rich Mg and 
Si isotopes using the HFB theory and exact angular momentum projection. 
In the calculations the finite range density dependent Gogny force has been 
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Table 2 

Transition probabilities B(E2,0 + — > 2 + ) in e 2 fm 4 between different configurations 
in 30,32,34 Msr and 32,34,36,38^ Ag in the prev j ous ta bi e q (P) stands for the oblate 

(prolate) configuration. The experimental data is taken from [4] (Si) and from [3] 
(Mg). 





0+ -» 2+ 


0+^2+ 


0+-2+ 


0+-2+ 


Exp. 


30 Mg 


182.11 


39.44 


0.88 


8.38 




32 Mg 


593.24 


15.53 


2.26 


6.52 


454±78 


u Mg 


549.21 


14.74 


2.66 


10.51 




32 Si 


402.09 


54.19 


4.26 


50.71 


113±33 


u Si 


227.18 


62.51 


6.81 


39.16 


85±33 


m Si 


232.22 


45.05 


6.43 


28.31 


193±59 


38 Si 


418.91 


35.18 


4.43 


69.80 


193±71 



Table 3 

Ground state energies in MeV as compared to the experimental results. The quan- 
tities SE I=0 = E Exp ~ E I=0 - E c .e. and 5E H fb = E Exp - E H fb - E c .e. are 
also presented. The quantity Ec.e. stands for the HFB Coulomb exchange energy 
computed in the Slater approximation. The experimental data are taken from [26]. 





Ehfb 


Ei=o 


Ec.e. 


EExp 


SE I=0 


SEhfb 


30 Mg 


-235.01 


-237.90 


-4.29 


-241.63 


0.56 


-2.33 


32 Mg 


-244.00 


-245.62 


-4.22 


-249.68 


0.16 


-1.46 


34 Mg 


-248.61 


-252.09 


-4.21 


-256.58 


-0.28 


-3.76 


32 Si 


-262.55 


-265.70 


-5.13 


-271.41 


-0.58 


-3.73 


u Si 


-276.32 


-277.59 


-5.06 


-283.42 


-0.77 


-2.04 


36 Si 


-284.15 


-287.22 


-5.04 


-292.01 


0.25 


-2.82 


38 Si 


-291.66 


-294.95 


-4.98 


-299.50 


0.43 


-2.86 



used. The results for the excitation energies + — 2 + and B(E2,0 + — > 2 + ) 
transition probabilities obtained from the angular momentum projected wave 
functions are in reasonable agreement with the experiment. The analysis of 
the projected energy surfaces and also the discrepancies found between theory 
and experiment indicate that configuration mixing is an important ingredient 
in these nuclei. Work is in progress in order to incorporate such configuration 
mixing. 
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the Spanish Instituto de Cooperacion Iberoamericana (ICI). This work has 
been supported in part by the DGICyT (Spain) under project PB97/0023. 
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